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On a Method for Determining the Non-Stationary 
State of Heat in an Ellipsoid. 



By Bibhutibhusan Datta. 



Introduction. 

1. The first writer, who attempted, with some success, the problem of the 
determination of the non-stationary state of heat in an ellipsoid with three 
unequal axes, was E. Mathieu,* who showed how the problem could be reduced 
to the solution of certain ordinary linear differential equations. But he found 
these equations to be so unmanageable that he contented himself with approx- 
imating to their solutions for the special case of an ellipsoid of revolution. 
Prof. C. Niven improved upon the results of Mathieu in certain respects in an 
interesting memoir,! entitled "On the Conduction of Heat in Ellipsoids of 
Bevolution." 

In the present paper, I propose (1) to obtain the chief results of Prof. 
Niven by an entirely different method, and (2), to show how this method can 
be applied to the case of the ellipsoid with three unequal axes, to obtain similar 
results which are believed to be new. It may be noted here that, in Art. 6, I 
point out a mistake in Prof. Niven's memoir. 

Preliminary Remarks and Definitions. 

2. Let the initial temperature of the ellipsoid be f(x, y, z) and let its 

x % u 2 z 2 
boundary, viz., —^-f — + -2=1, be maintained at temperature zero. Then 
a o c 

the required temperature V(x, y,z,t) is such that 

dV __ d 2 V d 2 V d 2 V 
dt ~ dx 2 + dy 2 + dz 2 ' { ' 

* Cours de physique mathematique, Ch. IX. t Phil. Trans., Vol. CLXXI (1880) . 
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the units being so chosen that the diffusivity is 1, 

V=0 on the boundary, (2) 

V=f(x,y,z) when t = 0. (3) 

Thus V can be expressed as a sum of terms of the form Ae~ XH W(x, y, z), 
where the normal function W satisfies the equation 

d 2 W d*W d 2 W 

&C 2 + dy* + dz* ~ ' 

and vanishes on the boundary; and the constants A are so chosen that the 
initial condition (3) is satisfied, so that f(x, y, z) ='LAW. 

When a = b=c so that the ellipsoid becomes a sphere of radius a, an 
appropriate normal function is 

S n (?i.r)PZ (cos 6) cos m<p, 
where 

S n {*>) = (-l) n ^2 X ~ iJ « + i(*)> p n(oos B) = sm m e-^, (i being cos 6, 

and X is a root of the equation 8 n {"Ka) =0. Throughout the present paper I 
will represent by W™ the normal function corresponding to the function 
S n (?»r)P™(cos 0) cos m$, and denote W° n by W n . 

I proceed now to obtain the functions W of various types. 

W For Ellipsoid of Revolution. 

3. Let e denote the eccentricity of the ellipsoid, then neglecting e 4 and 
higher powers, the equation of the ellipsoid can be written as 

r = ajl— £e 2 + ,}e 2 P 2 (cos0)f, i.e., r=a(l+eP 2 (cos 0) ), 
where a=a(l — J e 2 ) and g=^e 2 . Now assume that 

W = S {*r) +sXM t S t (lr)P t (cos $) , 

M t being an unknown constant to be determined. Then, evidently, W satisfies 
the partial differential equation (I) ; and to satisfy the boundary condition we 
must have 

O=S o (Xa)+£a^^P f +«2MA(;ia)P ( (coB0), 
since e 2 and higher powers are neglected. This equation must hold for all 
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values of cos 6. Therefore, equating to zero the coefficients of the various 

zonal harmonics, we get 

S {M) =0, (1) 

M 2 S 2 fr«)+ a d8o ^ a) =0, (2) 

and all the other M's are zero. Hence the required expression for W , in 
terms of a and e, is 

dS (Xa) 

(X 5 

W =S (lr)-ie* ™ S,ftr)P t (oosd), 

where % is given by the equation (1). But the general solution of the equation 

(1) is known to be 

3>a=t7t, 

i being any integer. Hence 

ha=in+$ ine 2 . 

4. In order to obtain a closer approximation to W , we will retain e 4 and 
neglect the sixth and higher powers. Thus the equation of the ellipsoid is 

r=a[(l-|e 2 -f 5 e 4 ) + (|e 2 +£e*)P 2 + 3 ^ 5 e 4 P 4 ], i.e., r=/3[l+<rP,+<rP 4 ], 

where /9=a(l-|e»-&e«), <x=|e»+SS«*» and T=£e«. 

Let us assume that 

dSpj^a) 

W =S (Xr)-ie 2 -^-^^S^r)P 2 (co S e)+^^N t S t (^r)P t (cos6), 

N t being an unknown constant to be determined. Expanding by Taylor's 
Theorem, we have 

8.{Xr) =8,m +3™^ (*P 2 +tP 4 ) + hP^JT 1 <<*.+***)'+ 
and S t (Xr) = S 2 &8) +8 d8 ^ 8) (ctP 2 +tP 4 ) + ...., 

when r=/?(l+crP 2 +TP 4 ). 

Again (P 2 ) 2 =i^+fP 2 +l. 

18 
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Hence we must have 

dS (Xa) dS (U) 

+TS2f,S,(V)P.- 

t=l 

This equation must hold for all values of cos 6. Therefore, equating to zero 
the coefficients of the various zonal harmonics, we get 

dS (*a) 

dS (Xa) 

■zN 2 b 2 {Z,p) + ap — g-g — ' + T a P — 9^2 ¥ e S % (*a) — **"***' 

a dS (Xa) 
2 e'er 8< * V^^ -Q (2) 

'tN i b i (AP) + tP gTj t-35 <T i y g^2 ■ 

'e'er * da -fl ^^ -0- (3) 

and all the other N's are zero. On substituting the values of 3, a and t in 
terms of a and e, the above equations take the forms 

3fl,(;ta) 
i 4 9a 3/S 2 (Xa) ... 

-- e -^^r a — a^-= ' (4) 

3£ (/la) 



3S JSI 2 b 2 (*a) + 21 a da - 21 o da , +63' £ i(Jla) a 9a ~ °> ( 5 > 

9^0 (Xa) 
^ A(M + ,ai^) +s> .~>_..^ r .a^)=0. ,6, 
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Hence the required expression for W is 

dS (Xa) 
W =S (lr) - | e 2 ^ 8 t (Xr)P t (G0B 6) + £ eW 2 S 2 (ar)P 2 (cos 0) 

+ &eWA(ar)Pi(eos0), 

the values of N 2 and N t being given by the equations (5) and (6), and the 
value of X being determined from the equation (4). 

5. I will now proceed with the solution of the equation (4) of the pre- 
ceding article. If we neglect e altogether, the equation reduces to S (Xa) =0, 
whose roots are given by %a=in. Therefore let the full value of /la be 

in+l 1 e l +l 2 e i + , 

l x , l 2 being unknown quantities which are to be determined. S n can be expanded 
in a series containing a finite number of terms so that 



o.ia»~-l 12 2V 



,n'(n'— l')(n'— 2')(n'— 3') 1 "I . / , nn\ 

, Tn' 1 n'(n'— l')(n'— 2') 1 , "I / , nn\ 

+ LT' 2^ O^ 2V + ••••J C0S ( a;+ ^/ , 

where n' stands for n(n-\-l) and (n'—f) for (%— *)(«+2+l). Hence we 
have the following: 

I e i II I 2 \ 
S (Xa)= -l-cosw+ t 1 -^ cos wi-e*+ , 

o /«> \ ^ • , 3# (/la) • h ' . , , 
c 2 (xa)=— -Ti— scos w-f- , a — ^ — -=cos*7i— t^ cos in- e 2 -f , 

a M^). a M^l^ 2(Aa) = | (iV -9)cosm+ 

Therefore, substituting the above expressions for S Q (Xa), etc., in the equation 
(4) of the preceding article, and equating to zero the coefficients of e 2 and e 4 , 
we obtain finally 

Xa-in-\- \ine* + ^ (iV+27)e 4 + 
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Comparison with the Results of Niven. 

6. By neglecting e 6 and higher powers, Prof. C. Niven has obtained an 
expression * for the parameter X, viz. : 

la=in+ |^e 2 + ^ (»V+27)e 4 + 

which differs from that obtained by me in the preceding article only in so far 
as the coefficient of e* is concerned. After carefully going through Prof. 
Niven's calculations, I find that the mistake of Prof. Niven must be attributed 
to some inadvertence on his part, as, by repeating his process, I get the correct 
result. 

From the identity of the values of X, it follows at once that my expression 
for W and Prof. Niven's must be identical. For, as is well-known, for the 
same value of A, there can not be two different solutions of the equation 

av a_¥ av 

both of which vanish on the boundary. 

W„ For Ellipsoid of Revolution, n>0. 

7. Let e 4 and higher powers be neglected, so that the equation of the 

ellipsoid is the same as in Art. 3, viz.: r=a\l-\-sP 2 (cos 0) \. Then assume 

that co 

W n =S n (lr)P n (cos 6) +eZ'H t S t (lr)P t (cos 6), 

00 

where the H's are unknown constants, to be determined, and S' refers to all 

(=0 

the values of t except t=n. Thus it is evident that W n satisfies the partial 
differential equation (I), and it remains to find the values of the constants H's 
so as to satisfy the boundary condition. Since e 2 and higher powers of e are 
to be neglected, we have, on putting r=a(l+eP 2 ), 

8 n ^r)=S n (Xa)+ea^^P 2i 
also, we have 

Pi ' Pn = #n+2Pn+2 + Cn-£ > » + -D„-2-F > n-2> 

where 






3 (n + l)(n + 2) n _ n(n + l) 

(2n+3)(2n+l) 

n(n — 1) 



2 (2n+3)(2n+l)' * (2n+3) (2»— 1) ' 



^»-2- 2 (2n + l)(2n— 1)' 



* hoc cit., p. 145. 
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Hence we must have 

0=S n ^a)P n +ea d8 ^ a) [B n+2 P n+i +C n P n +D n „ 2 P^l +ell'H t S t (^)P t . 
da, *=o 

This must be true for all values of cos 0. Therefore equating to zero the 
coefficients of the various zonal harmonics, we get 

S A Xa)+ea^^C n =0, (1) 

H n+2 S n+i (Ka) + a ^^B n+2 =0, (2) 

ff„_A- 2 (JW) + a^^D„_ 2 =0, (3) 

and all the other H's are zero. Thus the unknown constants are determined 
and the required expression for W „ , in terms of a and e, is 

a dSJMl 
W n =S n (*r)P- i e 2 8n J ( l a) S n+ >(*r)P n+2 B n+2 

a dsjMi 

~~ * ** S^(Xa) /Sf »-*( ar ) P «- 2l) »- 2 ' 
X being a root of the equation (1). 

But, in terms of a and e, this equation is 

6 » (aa) (2n+3)(2n-l) ea da ~°' (4) 

Therefore, if x be a root of the equation #„(#) =0, the corresponding solution 

of (4) is given by 

n*+n- 



Za=xll + 



n—l) J 



(2»+3) (2»- 

TF™ For Ellipsoid of Revolution. 
8. Let e 4 and higher powers be neglected, and assume that 

W%=S n (%r)P% (cos 6) cos m$+e 2' /f^O^Pf (cos 0) cos m$, 

t~m 

CO 

where 2' refers to all values of t from m up to infinity except the value t=n, 
and If is an unknown constant to be determined. Then it is evident that W™ 
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satisfies the partial differential equation (I). Now, putting r=a(l+eP 2 ), we 
get 

SAW =S n (la) +ea d8 ^ a) P 2 ; and P,-P7=B^P^ 2 +C^P:+D^P^, 

where 

7?m ~A (n—m+l)(n—m + 2) _ n(n+l)—3m* 

M+2-2 (2n+3)(2n+l) ' * (2»+3) (2n— 1) ' 

_ A (n+wp (n+<n— 1) 
B - 2 ~ 2 (2»+l)(2n— 1) ' 

Hence, from the boundary condition, we must have 

0=S n (te)PZ cos m»+6a a ^ a) [Br +1 P^ I +CrPr+D?_^L.] 



X cos m$+eI>'I?S t (%a)P? cos m$. 



t=m 



Therefore, equating to zero the coefficients of the various surface harmonics, 
we get 

sjk)+ w ^cr =o, 



da 

and all the other If's are zero. 

Thus the required expression for W% in terms of a and e is 



(1) 
(2) 
(3) 



w;= 



a- 



da 



S n {%r)P™(cosQ)-\e* 8 J S n+2 (kr)P^ 2 (eo3d)B^ 



a 



ie ! 



»+2(^a) 
da 



S B _ 2 (ar)Pr_ 2 (cos0)D:r_ 2 



cos w^), 



3 ^»-2(^) 

where X is a root of the equation (1). But, expressed in terms of a and e, 
this equation becomes 

s , (;u) _ (f+"-i)+"> ,,ag#gl=o. ( 4) 



Therefore, we obtain 



Aa 



-{ 



(2n+3)(2»— 1)"~ da 
i , (w 2 +n— l)+w z ^\ 



(2n+3)(2»— 1) 
corresponding to the root x of S n (x) =0 
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W For Ellipsoid with Three Unequal Axes. 
9. Let e x and e 2 be the eccentricities of the two principal diametral sec- 

SC 2 iP" £ 2 

tions of the ellipsoid — + 4« +-j =1, a>b>c, by planes passing through 

Qi C 

the major axis. Then, neglecting the fourth and higher powers of e x and e 2 , 
the equation to the ellipsoid is 

r=a[l- 1 (el+et) (P -P 2 ) + f 2 (e|-e?) Pi cos 2$>], 

i.e., r=y(l + e 1 P 2 +e 2 P 2 ! cos2$), 

where <y=a\l— | (ef+e|) \, 6 1 =^(e 2 1 +e 2 2 ), and e 2 =^(e|— e\). 

Now assume that 

W =S (Zr) + 2' 2 B mj< # { (/lr)Pr (cos 0) cos m$, 

oo 

where B m>t is an unknown constant to be determined and 2' refers to all values 

t=m 

of t from m up to oo, except £=0. Then it is evident that W satisfies the 
partial differential equation (I). To satisfy the boundary condition, we must 
have 

0=S (ty) +y d8o ^ y) (eiP 2 +e 2 Pi cos 2$) + 2' 2 R m , t 8 t (ty)P? (cos 0) cos mf, 

' ' Oy t=m m=0 

the unknown constants B's being assumed to be of the same order as ej. and e 2 . 

Therefore, equating to zero the coefficients of the various surface harmonics, 

we get 

S (Xy) =0, (1) 

Bo, i S i ^y)+e 1 y^^-=0, (2) 

R>,>8 2 (ly)+e 2 y^2±=0; (3) 

and all the other R's are zero. Thus the R's are determined, and we get 
finally the required expression for W Q to be 

dS (?ia) 



W Q =S (%r)-± (el+et) g ,,„, S 2 (lr)P 2 (cos 6) 



a 



da 
8 t (Xa) 
dS (ha) 



+ h ( e i~ e i) g^j 8 2 &r)P\ (cos 6) cos 2*, 
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where X is a root of the equation (1). But the roots of the equation (1) are 

given by 

ty=irt. 

Hence, in terms of a and the eccentricities, we have 

to=in+ ■§■ (ef+e|)«i. 

Conclusion. 

10. I conclude this paper by pointing out that the results or Arts. 3-9 
admit of extension and generalization in various directions. For example, a 
procedure similar to that of Art. 9 will give us W? for the ellipsoid with three 
unequal axes. Also, denoting by W JT, the normal function corresponding to 

S n (%r)P% (cos 0) sin m$, 

it is obvious that, for the ellipsoid of revolution as well as for the ellipsoid 
with three unequal axes, W% can be obtained in the same way as W?. 

Univebsity of Calcutta, 1917. 



